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duct cross section in the yz -plane changes from circular 
to nearly rectangular in the region 1.0 < xf D\ < 2.5. 
The duct aspect ratio (the ratio of the major to minor 
axis lengths at xjD\ > 2.5) is 3.0. The cross-sectional 
areas at x/D\ <1.0 and xjD\ > 2.5 are equal. In 
the region 1.0 < xjD\ < 2.5 the cross-sectional area 
increases as much as 1.15 times the area at xj D\ < 1.0. 

Experiment 

The inlet diameter of both the S-duct and tran- 
sition duct was 10.24 cm. The S-duct is larger and 
geometrically similar to the duct studied by Vakili el 
al. 2 The transition duct is identical to that of Davis 
and Gessner. 4 The inlet conditions for both ducts are 
summarized in Table 1. For the transition duct ex- 
periment with swirling flow, the maximum swirl an- 
gle was 15.6°. A swirl generator was used to produce 
solid body rotation with minimal downstream distur- 
bances. Additional information about the experimental 
pro gram is contained in Refs. 7, 8. 

Detailed measurements of velocity, total pres- 
sure, and static pressure were made at cross -stream 
planes inside both ducts with calibrated three- and 
five-hole probes. Data were acquired at about 530 
locations near the S-duct inlet (at $jD\ = —0.5) 
and at nearly 1220 locations near the S-duct exit (at 
s/D\ = 5.73). In the transition duct, measure- 
ments were made in four cross -stream planes located 
at x/D\ = 1.49, 1.99, 2.55, and 3.93. The measure- 
ments were acquired at approximately 480 locations in 
each measurement plane. Surface oil film visualization 
and surface static pressure measurements were also ob- 
tained for both ducts. 

Table 1 Experimental Inlet conditions. 


S-duct Transition duct 
No swirl Swirl 


-^centerline 

0.60 

0.50 

0.35 

Re Dl x 10 -6 

2.20 

2.09 

1.37 

$o. 95/^1 K 100 

3.04 

3.76 

10.24 


Computation 

The PARC3D computer program solves the full, 
three-dimensional Reynolds- averaged Navier-Stokes 

equations in strong conservation form with the Beam 
and Warming approximate factorization algorithm. The 
implicit scheme uses central differencing for a curvilin- 
ear set of coordinates. Ref. 9 describes in detail the 
theory and features of the PARC code. 

A Baldwin -Lorn ax 10 algebraic turbulence model 
was used. The model was modified to use only vorticity 
in the local boundary layer to avoid secondary flow ef- 
fects on the eddy viscosity. 11 The turbulence model was 
also modified for the swirling flow transition duct com- 
putation to remove the constant level of axial vorticity 
from the vorticity level in the boundary layer. For the 
S-duct computation a low- Reynolds number k — c tur- 
bulence model of Speziale et al. n was also used. The 
numerical solution algorithm is that of Nichols. 13 



s/Dj 

Fig. 3 Surface static pressure coefficient 
for S-duct flow, (a) k — e turbulence 
model, (b) algebraic turbulence model. 

The S-duct grid was composed of three blocks with 
grid distributions of 32 x 71 x 53, 69 x 71 x 53, and 
32 x 71 x 53. An H-grid of 129 x 11 x 15 was used 
in the center. Gridgen 14 was used to define the interior 
grid. The first point off the wall had an average y + of 
less than 1. The boundary conditions were no slip at 
the walls, total pressure and temperature specified at the 
entry plane, static pressure specified at the exit plane, 
and symmetry about the x z -plane. 

For nonswirling transition duct flow the grid mod- 
eled one duct quadrant and contained 97 x 51 x 53 
grid points. The inlet total temperature and pressure 
were specified locally. At the duct exit, a constant static 
pressure was specified and density and velocity were 
extrapolated from the interior. The downstream static 
pressure was adjusted to match the measured surface 
static pressures at the inlet. The swirling flow transition 
duct computation was performed with a 97 x 51 x 97 
grid that modeled two duct quadrants. New boundary 
conditions were incorporated in the PARC code to per- 
mit swirling inlet flow and to accommodate the 180° 
rotational symmetry of the flow. For the swirl case, the 
inlet boundary condition used the experimentally mea- 
sured inlet flow angles. 



s/Di 

Fig. 4 Skin friction coefficient for 
S-duct flow, (a) k — e turbulence 
model, (b) algebraic turbulence model. 
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(a) Experiment Computation 


(b) Experiment Computation 

Fig. 5 Exit plane total pressure coefficient 
for S-duct flow, (a) k — e turbulence 
model, (b) algebraic turbulence model. 

Results and Discussion 

The total and static pressures in this section are 
presented as pressure coefficients given by Eq. (1). 
The pressures p 0 and p represent the local total and 
static pressure. The reference variables, subscripted cl 
(centerline) or wall, were evaluated at one-half inlet 
diameter upstream of the S-duct or transition duct inlet 
(at s/Di = -0.5 or z/Di = -0.5). 

Cpo = , Cp= ( 1 ) 

PO,cl Pwall P0,cl ~ Pwall 

Diffusing S-Duct Results 

The computed and experimental surface static pres- 
sure distributions are shown in Fig. 3. Both the alge- 
braic and k — c turbulence model results agree with the 
experimental values upstream of the experimentally de- 
termined separation (shaded region in Figs. 3 and 4). In 
the separation zone, the k — f model agrees better with 
experimental values, and downstream of separation both 
models begin to agree with the experimental data. Pre- 
dicted and experimental skin friction coefficient values 
are plotted in Fig. 4. The agreement is considered rea- 
sonable. The experimentally measured region of flow 



Fig. 6 Law of the wall velocity profiles for 
flow near S-duct exit, (a) k — c turbulence 
model, (b) algebraic turbulence model. 

separation is 2.02 < sjD\ < 4.13. The k - e and the 
algebraic turbulence models predict separation between 
2.59 < s/Dx < 4.25 and 2.69 < s/Di < 4.25, re- 
spectively. The computational values were determined 



(a) Experiment Computation 



(b) Experiment Computation 


Fig. 7 Exit plane transverse Mach vector 
components for S-duct flow, (a) k — c turbulence 
model, (b) algebraic turbulence model. 
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Fig. 8 Near surface flow visualization, 
(a) k — e turbulence model, (b) algebraic 
turbulence model, (c) experiment. 


from locations of negative axial velocity at the grid point 
nearest the wall along (j> — 180°. 

Computational and experimental total pressure con- 
tours near the S-duct exit (sf D\ = 5.73) are compared 
in Fig. 5. The predicted region of diminished total 
pressure is smaller than the measured region for both 
turbulence models. This discrepancy is possibly due to 
a turbulence modeling deficiency to account for with 
strong cross flow (three-dimensional) effects or a mod- 
eling deficiency of artificial viscosity in the boundary 
layer separation region. The grid resolution should be 
adequate as the distance of the first grid point off the 
wall is approximately y + of 1- 

The predicted velocity profiles for both turbulence 
models near the S-duct exit arc shown in wall coordi- 
nates in Fig. 6. Near the S-duct exit the vortices have 
convected low velocity fluid away from the wall and the 
boundary layer in this region, particularly at <j> = 170°, 
departs significantly from this law. 

Experimental and computed transverse velocity 
vectors are shown in Fig. 7 near the duct exit. The re- 
sults show the presence of large counter-rotating vortices 
in the lower half of the cross section. Both computations 
are in qualitative agreement with the experiment. 

The predicted streamlines in the region of flow sep- 
aration are compared with the surface oil film visualiza- 
tion in Fig. 8. The predicted flow fields are in general 
agreement with the data, and the k — e turbulence model 
streamlines better compare to the experimental stream- 
lines. 


Circular-to-Rectangular Transition Duct Results 

The numerical and experimental values of surface 
static pressure for both the nons willing and swirling flow 
cases are plotted in Fig. 9. The prediction is very good 
indicating that the aerodynamic blockage is correctly 
predicted. Experimental and computed contours of the 
total pressure coefficient at x/ D\ — 2.55 are shown in 
Fig. 10 (a) for flow without inlet swirl. The concentric 
contours extending inward from the duct side walls are 
regions of diminished total pressure that result from 



Fig. 9 Surface static pressure coefficient 
for transition duct flow, (a) without 
inlet swirl, (b) with inlet swirl. 


side wall vortices convecting low total pressure fluid 
away from the duct surface. This is the most upstream 
location where the effects of the side wall vortices are 
noticeable in the experimental data. The numerical 
results predict the same flow structure at x/ D\ = 2.55, 
but the region of diminished total pressure is not as large 
as the experimental results reveal. 

Total pressure contours at x/Z>i =3.93, Fig. 10 
(b), also show the predicted region of diminished total 
pressure does not extend as far from the side walls as 
was experimentally observed. Comparing experimen- 
tal and numerical surface flow visualization indicates 
that the point of formation of the vortices is correctly 
predicted. This suggests that the turbulence model, re- 
stricted to wall bounded shear flow, can not account for 
turbulent mixing in the vortex region. The agreement 
outside the region affected by the side wall vortices at 
xjD\ — 2.55 and 3.93 is excellent. 




Computation 



Computation 


Fig. 10 Total pressure coefficient for 
transition duct flow without inlet swirl, 

(a) at x/D l = 2.55, (b) at xjD x = 3.93. 
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Fig. 11 Total pressure coefficient for 
transition duct flow with inlet swirl, (a) 

at x/D x = 2.55, (b) at x/D x = 3.93. 

The total pressure at x/D x = 2.55 and 3.93 for 
flow with inlet swirl are shown in Figs. 1 1 (a) and (b). 
With swirl the flow is symmetric with respect to 180° 
rotations about the x-axis. Regions of diminished total 
pressure near the duct side walls (that were produced 
in the nonswirling flow by the side wall vortices) are 
absent. However, surface oil film visualization, Ref. 7, 
indicates cross flow near the duct corners in the down- 
stream region of the duct. The effect of the cross flow 
on the total pressure coefficients is visible, particularly 
in the upper left (experiment) and lower right (com- 
putation) quadrant at x/D\ = 2.55. In general, the 
agreement between the experiment and computation is 
better for the swirling case (no side wall vortices) than 
the nonswirling case (side wall vortices). This is most 
apparent at x/Di =3.93 where the comparison in Fig. 
11 (b) is noticeably better than Fig. 10 (b). 

Conclusion 

The PARC3D Navier-Stokes code has been used to 
compute flow through a diffusing S-duct and a circular- 
to-rectangular transition duct with and without inlet 
swirl. The objective of these calculations was to assess 
the ability of the PARC code to accurately predict flows 
through propulsion system ducts. The S-duct computed 
flow fields are generally in good agreement with the ex- 
perimental data. However, both turbulence models un- 
derpredict the length and angular extent of the boundary 
layer separation, and the predicted separation occurs one 
half sj D\ late. Neither algebraic nor k — e turbulence 
model adequately accounts for strong secondary flows 
with separation. The transition duct flow field matched 
the experimental total and surface static pressure coeffi- 
cients well. The agreement appeared better for flow with 
inlet swirl, where the pairs of counter-rotating vortices at 
the duct exit were absent. For attached flow, PARC3D 
has demonstrated reasonable accuracy. For flows with 
strong cross flow and or with separated boundary layers, 
the modeled turbulence and or artificial viscosity should 
be improved. 
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